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THE EFFECT OF WALL FRICTION ON MAGNETOHYDRODYNAMIC 
GENERATOR PERFORMANCE 
by Allan R. Bishop 
Lewis Research Center 

SUMMARY 

The effect of wall friction on magnetohydrodynamic generator performance is deter- 
mined by introduction of a wall friction factor into the one-dimensional generator equa- 
tions. This addition should be useful in improving generator analysis and determining 
optimum generator geometry. The curves presented in this report can be used to deter- 
mine the effects of changes in wall friction and generator geometry on generator per- 
formance. Wall friction has an increasing effect as the Mach number increases and a 
decreasing effect as the pressure drop across the generator increases. 


INTRODUCTION 

Magnetohydrodynamic (MHD) power -generator concepts can be divided into two 
principal classes: (1) large, high-power, ground-based generators, and (2) lower -power 
generators for space power applications. In ground-based generator concepts, the chan- 
nel volumes are large, and magnetic field strengths are high. Wall friction effects on 
generator efficiency are therefore small, and the correlation between theory and experi- 
ment is generally good (refs. 1 and 2). The weight of the MHD generator and its magnet 
are not serious problems in ground-based generators. 

For space power applications, the weight of generator and magnet can be important. 
Hence, a compact, lightweight system is needed, which implies small channel cross 
sections where friction effects can be important. Until now, most MHD generator anal- 
yses (e. g. , ref. 3) for space power applications have neglected the effects of wall fric- 
tion. 

Nonuniformities in the fluid flow field and in electrical current distribution can also 
affect MHD generator performance. The magnitude of these effects may be smaller or 
greater than the effects erf wall friction, depending on the generator configuration. No 
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attempt has been made to include nonuniformities in this analysis. The generator design, 
basic assumption, and solution techniques used in this report are similar to those of 
reference 3, except that a wall friction factor has been introduced into the one- 
dimensional momentum equation. 


ANALYSIS 

The generator is assumed to be of the Faraday segmented type, with infinitely fine 
segments so that no current flows in the Hall direction. The working fluid is assumed 
to be a perfect gas, with zero viscosity and thermal conductivity. A friction factor is 
introduced to correct for the drag of the generator walls (ref. 4). 

The following procedure is used in this analysis: 

(1) The one -dimensional flow equations are written in nondimensional form. 

(2) The assumptions of uniform power generation and uniform Joule dissipation along 
the generator are used to reduce the number of unknown variables. 

(3) The variation of the cross-sectional area is determined by requiring minimum 
generator volume for a given total pressure ratio. A variational method with Lagrange 
multipliers is used to minimize the volume. 

(4) With the area variation known, the generator performance can be found as a func- 
tion of inlet and load conditions. 

All symbols are defined in the appendix, and all quantities are expressed in the 
International System of Units. 

The five one -dimensional equations for the flow are 

Continuity: 


— (pua) = 0 (1) 

dx 


Momentum: 


du , dp f F pu 

pu + — Hr jB + - 

dx dx 2r^ 


= 0 


( 2 ) 


Energy: 


pu — + pu 
dx 


2 - je = 0 

dx 


( 3 ) 
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where 


Ohm's law: 




y * 1 P 


(4) 


j = a(uB + E) (5) 

The last term in equation (2) is the correction for wall friction losses, and is the only 
addition to the flow equations in reference 3. 

These equations can be written in nondimensional form by using the following defini- 
tions; 


A = 


a 




p l 


X * 


p 1 U 1 


u = — 

Ui 




B 


P = 


^l u l 


J = 


OjUiBj 


where the subscript 1 indicates channel entrance conditions. Equations (1) to (5) re- 
duce to 


gg + Agg + JAC* (fPiUi C = 0 


dX dX 




r 


- — (AUP) + U— + JAC [U - = 0 

1 dX dX V SC/ 


Define a new parameter F Q such that 


F o = 


f F PjUi 




( 6 ) 

(7) 



Note that r h is the hydraulic radius, which is one-fourth of the hydraulic diameter 
and is assumed to be constant. Equation (6) then becomes 

— + A — - + JAC + FU = 0 (8) 

dX dX ° 


It is now assumed that 


JC = J 


ent 


and 


J 

SC 


J 


ent 


where J gnt is a constant. The use of these assumptions results in uniform power 
generation and uniform Joule dissipation along the generator. 

The constant ^ent is related to the load parameter K by the equation 


where K is the ratio of actual output voltage to open circuit voltage. The value of J ^ 
is 0 at open circuit and 1 at short circuit. 

With these two additional assumptions, equations (7) and (8) reduce to 


Y 


y 


_ J. (AUP) + U “ + J 
1 dX dX 


ent 


A(U - J ent ) - 0 


(9a) 


and 


— + A — + J .A + F U = 0 (9b) 

dX dX ent 0 


The total pressure in dimensionless variables is 


where 



Let 

q T = (P X ) M 

and 


q * (P) M 

Then, solving for U and substituting into equation (9b) yields 

dU + A dP + jj ent + l F Q (q T - q)q l/M_1 jA dX = 0 (11) 

It is still necessary to specify the variation in cross-sectional area as a function of 
the distance along the generator length. It is important to maintain a high generator ef- 
ficiency for a specified total pressure ratio. The only energy loss considered in this 
analysis is Joule heating, which varies with generator volume and is assumed uniform 
along the duct. It seems reasonable, therefore, to minimize the generator volume for a 
specified total pressure ratio, since this will minimize the energy loss. The procedure 
has also been followed in reference 3. 

The volume r of the generator is given by the integral of the area 

f l 

t = ^ a dx 


Define a nondimen si onal volume V as 


V = 


*Wi 


-r 


A dX 


( 12 ) 
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Substituting equation (11) into (12) yields 


V = - 



dU + A dP 

- f 0 (q t - - q)q 1/M-1 


Define 


*9- +A -!® 


F = 


dq„ 


dq„ 




then V is 


V = - 


/ 


exit 


F dq n 


This is the integral to be minimized, subject to the constraints of equations (9) 
and (11). Eliminating A dX from equations (9) and (11) yields 


fdu + u_Sa 

\ 2 q T 


^ent< D - J ent> = [ J ent + ~ <«T ' OX^"' 1 ] 

r~ 

Uq T dq Uq dq T \ 

- q/ 


J L- L 

2(q T - q) \ 


2q dU + 


dU 


q T - q q x 

At the minimum value of V, the variation in V must vanish (ref. 5, p. 74). 

6q T 


6V = - — 3F 

«u-ii£ 

6q - - F 

2 3U' 

L 2a< >' 

L 2 




where 


dq T 

dq T 

The third term on the right side of equation (16) is zero, since the exit total pres- 
sure q T is fixed. The first two terms may be reduced to the constraint equation (15) 
with dq T set equal to zero, and their sum must be zero. Therefore, the integral term 
must also vanish. The minimization of the volume may then be handled as a fixed end- 
point problem with a constraint. With the introduction of the new variable 


v = 


U 


qq* “ q 


(17) 


and its derivative 


v* - d Y - U 1 

dq T q T - q 

the integrand of equation (14) becomes 


v 

q T - q 


(1 - q’) 


F = 


v’(q*r - q) - q’ - + v 

1 2 


'ent 


2F, 


4 


(18) 


(q T - q)q 


i/ (i-i 


and the constraint equation is 

f = [v'(q T - q) - ^ q’J ^ent (J ent + «*> F o v M% ‘ 

+ v [ J ent( J ent + f ) + ~ (q T ‘ q)(2<1 T ' q)q 1/M ^ ° 

Applying the technique of Lagrange multipliers (ref. 5, p. 129), introduce 
H = F + \(p, where X is the undetermined Lagrange multiplier. The integral 


(19) 


I will be 
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minimized if the equations 


3H __ _d_£H = 0 
3q dq T 9q' 

3H _ _d_ SH _ 0 l 
3v dq T 3v' 

(p = 0 _> 


( 20 ) 


are satisfied. Substituting in H and taking the appropriate derivatives gives three 
equations for q, v, and X: 



where 


a 


J ent . 2F o<«T ' «> 

2 2 

M 

jq x U - m) - qj 

H)] 

q ‘/M-2 


j ent+ ^(q T -q)q 1/, ‘' 1 

2 


( 21 ) 


+ XJ 


ent 


v(q T - q) 


ent 


2XF c vq T 




(q x - q)[q T U - m) - q]q 


l/p-2 


» - 1 [ J ent< J ent * « v » + 

b ' f* 

* 


y — " (^t * *l) 


J ent< J ent * <l v > + 
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vF ?t q T^ 1 " ^ ~ q ^ 2 " M)]q 1//M 


J ent + ^ q T “ q ^ ql//,X 


-2 

r 


Xv 2 <N + F 

2 M 2 


-2 [ 2 q|d - M) “ QQ t (3 - 2 4 ) + q 2 ]q 1/#i “ 2 


B* = 


4 


2F 0 (q T - Ih 1 ^- 1 ,^1/n-l 

lq T 


J ent + ~ ’ q)ql/4 “ 1 
® 4 

F v 


C* = - v 


T ent^ + J ent) + (q T " q)(2q T " q)q 


1 / 4-1 


The matrix equation (21) is singular and unless 


G = oC* + @B* - yA* = 0 


( 22 ) 


there are no solutions. Equation (22), which is required for the matrix to have a solution 
defines the Lagrange multiplier X. 

With the conditions G ~ 0 and dG = 0, the matrix can be rewritten as 



(23) 


The solution to this system of equations determines the variation of v, q, and X 
which minimizes the volume. Substitution into equation (10) gives the proper area varia- 
tion, and the volume can be determined from equation (12). 



RESULTS AND DISCUSSION 


Matrix inversion by the Gauss- Jordon reduction method was used to solve equa- 
tion (23). Typical values for the dimensionless volume and area ratio are shown as func- 
tions of total pressure ratio in figure 1, for an entrance Mach number of 1. 0 and J gn ^ 
of 0. 25. The solid lines are the solution for zero friction, F q = 0. 0, while the broken 
lines are the solutions for F Q = 0. 1. 

The curves for F Q = 0. 0 are identical to the curves in figure 7 of reference 3. The 
value F Q = 0. 1 was determined from the following typical entrance conditions: u^ = 800 
m/sec, pj = 0. 25 kg/m^, f F = 2. 5x10"^, r h = 2. 5 cm, = 1.0 tesla, = ohm/m. 

Since the area and volume are variables, a change in the friction can alter the area and 
volume for a given pressure drop. Note that the area increases slightly but the volume 
decreases as F Q is increased at constant pressure. 

In figure 2, the nondimensional generator length is plotted as a function of the 
total pressure ratio for two values oi F q , As might be expected, the same pressure 



Total pressure ratio, Py 

figure 1. - Dimensionless volume and area ratio as 
functions of total pressure ratio* Entrance Mach 
numlter, LOj dimensionless load current, 0,25, 
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0 A A .6 A 10 


Total pressure ratio, Py 

Figure 2. - Dimensionless generator length as a function 
of Iota! pressure ratio. Entrance Mach number, LG; 
dimensionless load current, 0.25. 

drop occurs in a shorter length as the friction forces are increased (F increased)* The 

increase in cross-sectional area (fig* 1) is not enough to compensate for the increased 

friction* As a result, the volume is smaller at the increased friction* 

In most MHD generators, an important design factor is the generator efficiency* 

Here the efficiency means the isentropic efficiency, the calculated change in enthalpy 

divided by the isentropic change in enthalpy between the same total pressure limits* 

Sutton and Sherman (ref* 6) discuss the effects of wall friction on generator efficiency 

in terms of a polytropic efficiency . The appropriate expression for rf given by 

P P 

equation (14. 98) in reference 6 is 
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where the notation of reference 6 has been changed to correspond to the notation used in 
this report. Since U, A, and M may vary as the pressure ratio changes along the 
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length of the generator, r? p is actually a local polytropic efficiency. The generator effi- 
ciency ri„, from reference 6, is 

O 




g 


1 - p 


n p (r-i)/y 

T 


i - p 


(y-i )/y 

T 


This efficiency depends on the total pressure ratio, both explicitly and implicitly, 
since ?j p also can depend on pressure ratio. In order to calculate the generator effi- 
ciency as a function of pressure ratio, it is necessary to include the polytropic efficiency 
variation. The effect of the variation in polytropic efficiency can be seen in figure 3(a). 
The generator efficiency is calculated for three cases: (a) a constant polytropic effi- 
ciency evaluated at the entrance condition, (b) a constant polytropic efficiency evaluated 
at the exit condition, and (c) the correct local polytropic efficiency evaluated at each 
pressure from the entrance to the exit. Curve c fits between curves a and b, just as 
one would expect. 



Total pressure ratio, Pj 
(a) Friction factor 0, 1. 



♦ 24 


Figure 3. - Generator efficiency as a function of total pressure 
ratio and friction factor. Entrance Mach number, 10; di- 
me ntlon less load current, 0.25. 
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The overall generator efficiency is therefore dependent on the pressure ratio, which 
in turn is dependent on generator length. The efficiency is also a function of hydraulic 
radius. The efficiency is plotted in figure 3(b) as a function of F Q (which is inversely 
proportional to the hydraulic radius). The curves labeled a, b, and c correspond to 
those in figure 3(a) . Again, with the use of the polytropic efficiency at each point, 
curve c fits between curves a and b, which are based on entrance and exit conditions. 

Thus, the efficiency is affected by both hydraulic radius and length, but a new non- 
dimensional number that combines both factors can be defined. Multiplying F Q by the 
nondimens ional generator length yields 


„ Vl»l 

j* 4 A T = — * — — * — 

*vA p » u * 


2r h 

where L is the generator length, fp is the friction coefficient, and r^ is the hydraulic 
radius. This parameter is plotted in figure 4 as a function of J for constant genera- 
tor efficiency. Curves are plotted for two pressure ratios, 0. 5 and 0. 1. At pressure 
ratios lower than 0. 1, the area ratio increases so rapidly (fig. 1) that a one -dimensional 
analysis no longer applies. In figure 4(a) the entrance Mach number is 0. 5, in fig- 
ure 4(b) it is 1. 0, and in figure 4(c) it is 2. 0. 

From the curves in figure 4, we can determine the effect of friction on the generator 
efficiency. Clearly, as the friction coefficient increases (or the generator length L in- 
creases, or the generator hydraulic radius r^ decreases), the generator efficiency goes 
down at constant J gn p To maintain constant efficiency, J gnt could be lowered, but with 
the penalty of decreased output power or increased generator size. 

As an example, the 10 -MW cesium -seeded generator mentioned in reference 7 has 
an inlet gas temperature of 2500 K at 10 atmospheres and is 1/2 meter long. Assuming 
a square cross section, an aspect ratio (length divided by square root of entrance area) 
of 3, an 80-percent efficiency, and operation near Mach 1, the required entrance mag- 
netic field for a frictionless generator is 8. 4 tesla. If the friction coefficient fp is in- 
creased to 0. 003, the magnetic field required to maintain constant efficiency increases 
to 10. 0 tesla. If the magnetic field is not increased, the generator efficiency will be re- 
duced to approximately 70 percent. A friction coefficient of 0. 003 is within the current 
state of the art, but some care is required to achieve it. 

Generator efficiency is much more sensitive to wall friction as the Mach number in- 
creases, as can readily be seen by comparing figure 4(a) with figures 4(b) and (c). To 
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maintain constant efficiency, a larger change in J gnt is required at Mach 2.0 than at 
Mach 0. 5 (fig. 5), for the same change in friction coefficient. The results of this anal- 
ysis indicate that smooth walls and proper geometry are much more important in a 
supersonic generator than in a subsonic one. 

The expansions to a pressure ratio of 0. 1 seem to be less sensitive to friction than 
the expansion to a pressure ratio of 0. 5. Note that in figure 4 the 0. 1 expansion has a 
steeper slope than the expansion to 0. 5. The increased area at lower pressures (fig. 1) 
and the resulting reduction in wall friction effects may explain this difference. 

None of the curves in figures 4 or 5 have been extended to J en £ = 0* because a solu- 
tion could not be obtained in this region. The mathematical difficulty was a zero value 
for the determinant in the matrix inversion procedure. Whether some physical phenom- 
enon was occurring or whether a mathematical instability of the analytical procedure was 
at fault could not be determined. For some, but not all, cases the problem was asso- 



Figure 5. - Dimensionless friction coefficient 
as a function of entrance current density 
ratio, with Mach number as a parameter* 
Efficiency, 0.75; total pressure ratio, 0,5. 
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ciated with an adverse pressure gradient in the flow. A similar situation was noted in 
the analysis of reference 1, but since the analytical methods are quite different, a di- 
rect comparison was not possible. 


CONCLUSIONS 

The addition of wall friction to the MHD generator equations should be useful in im- 
proving generator analysis. The curves presented in this report can be used to deter- 
mine optimum generator geometry and the effect of changes in wall friction and gener- 
ator geometry on generator performance. 

The wall friction seems to have a greater effect as the Mach number increases, and 
a smaller effect at large pressure drops across the generator. The latter effect may be 
the result of the area variation used in this analysis. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 30, 1972, 

112-02. 
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APPENDIX - SYMBOLS 


A 

A*, B*, C* 

a 

B 

C 

E 

F 



j 

K 

L 

P 

P 

q 

R 

r h 

u 

u 

V 

v 

X 

X L 


generator area ratio, a/a^ 

terms defined in eq. (21) 

generator cross-sectional area 

magnetic field 

magnetic field ratio, B/B ^ 

generator electric field 

integrand in minimum volume integral 

nondimensional friction factor 

wall friction coefficient 

parameter defined by eq. (22) 

function defined by eq. (20) 

enthalpy 

dimensionless current density 
dimensionless load current 
current density 
load parameter 
generator length 
dimensionless pressure 
pressure 

variables defined in eq. (10) 
density ratio p/pj 

hydraulic radius; ratio of cross-sectional area to perimeter 

velocity ratio u/uj 

velocity 

dimensionless generator volume 
variable defined in eq. (17) 
dimensionless length 
dimensionless generator length 



a, P,y terms defined in eq. (21) 
y ratio of specific heats 

6 variational symbol 

rj generator efficiency 

D 

77 polytropic efficient 

F 

X Lagrange multiplier 

j u function of specific heat ratio defined in eq. (10) 

p working fluid density 

2 electrical conductivity ratio cr/tr ^ 

a electrical conductivity 

r generator volume 

<p constraint equation defined in eq. ( 19 ) 

Subscripts: 

exit exit conditions 

L generator exit value 

T stagnation property 

1 entrance conditions 
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